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Abstract. The 11 -year solar cycle in solar spectral irradi- 35 
ance (SSI) inferred from measurements by the SOIar Ra- 
diation & Climate Experiment (SORCE) suggests a much 
larger variation in the ultraviolet than previously accepted. 

5 We present middle atmosphere ozone and temperature re- 
sponses to the solar cycles in SORCE SSI and the ubiquitous 
Naval Research Laboratory (NRL) SSI reconstruction us- 
ing the Goddard Earth Observing System chemistry-climate 40 
model (GEOS CCM). The results are largely consistent with 
10 other recent modeling studies. The modeled ozone response 
is positive throughout the stratosphere and lower mesosphere 
using the NRL SSI, while the SORCE SSI produces a re- 
sponse that is larger in the lower stratosphere but out of phase 45 
with respect to total solar irradiance above 45 km. The mod- 
15 eled responses in total ozone are similar to those derived 
from satellite and ground-based measurements, 3-6 Dob- 
son Units per 100 units of 10.7-cm radio flux (F 10 . 7 ) in the 
tropics. The peak zonal mean tropical temperature response 50 
using the SORCE SSI is nearly 2K per 100 units F 10 . 7 — 
so 3 times larger than the simulation using the NRL SSI. The 
GEOS CCM and the Goddard Space Flight Center (GSFC) 
2-D coupled model are used to examine how the SSI so- 
lar cycle affects the atmosphere through direct solar heating 56 
and photolysis processes individually. Middle atmosphere 
25 ozone is affected almost entirely through photolysis, whereas 
the solar cycle in temperature is caused both through direct 
heating and photolysis feedbacks, processes that are mostly 
linearly separable. Further, the net ozone response results e0 
from the balance of ozone production at wavelengths less 
30 than 242 11 m and destruction at longer wavelengths, coinci- 
dentally corresponding to the wavelength regimes of the SO- 
Lar STellar Irradiance Comparison Experiment (SOLSTICE) 
and Spectral Irradiance Monitor (SIM) on SORCE, respec- 65 
tively. A higher wavelength-resolution analysis of the spec- 
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tral response could allow for a better prediction of the atmo- 
spheric response to arbitrary SSI variations. 


1 Introduction 

Recent measurements of the solar spectral irradiance (SSI) 
made by the the SOIar Radiation & Climate Experiment 
(SORCE) mission (Rottman, 2005) have brought into ques- 
tion long-held views of how the SSI varies as a function of 
(he 1 1-year solar cycle. Based on SORCE observations dur- 
ing the latter part of the declining phase of solar cycle 23 
(2004—2007), Harder et al. (2009) suggested that the mag- 
nitude of the solar cycle variation in the ultraviolet (UV), 
from 200 to 400 nm, is several times larger than previously 
thought. Further, they posited that the SORCE SSI variation 
in the visible and near-infrared (NIR) is out of phase with the 
UV and total solar irradiance (TSI) (i.e., the visible and NIR 
flux actually increased approaching the last solar minimum). 
The larger UV variation has significant implications for the 
response of the stratosphere and mesosphere and potentially 
climate to the solar cycle. 

Several studies (Cahalan et al., 2010; Haigh et al., 2010; 
Merkel et al., 2011; Ineson et al., 201 1) suggest that atmo- 
spheric models forced with the SORCE SSI lead to better 
agreement with the solar cycle response inferred from obser- 
vations than simulations based on the widely accepted Naval 
Research Laboratory (NRL) SSI reconstruction derived from 
long-term solar observations (Lean, 2000), implying that the 
solar cycle in SSI derived from SORCE data is more consis- 
tent with observations. Garcia (2010) and Matthes (201 1), 
however, have noted caveats with both the completeness of 
the model calculations and the short data records used. In ad- 
dition, Morrill et al. (2012) and DeLand and Cebula (2012) 
have more directly examined the SSI measurements them- 
selves in comparison to other coincident and historical obser- 
vations. Although the SORCE SSI measurements have been 
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found to be internally consistent and agree with other ob- 
servations on short time scales (e.g., 27-day solar rotation), 
these studies question the veracity of the 1 1 -year solar cycle- 
derived trends. 125 

In this paper we present results of simulations made with 
two state-of-the-art chemistry-climate models (CCMs), the 
Goddard Earth Observing System (GEOS) 3-D CCM and the 
Goddard Space Flight Center (GSFC) 2-D model, using both 
the NRL and SORCE SSI. We examine the individual mech- 130 
anisms of atmospheric response to solar cycle variations in 
SSI through direct solar heating and photolysis. We also ex- 
plore how middle atmosphere ozone responds to solar cycle 
variations in different parts of the SSI spectrum and what im- 
plications can be drawn regarding the inferred SORCE solar 135 
cycle in SSI. 


2 Background 

140 

Temporal variations in the solar energy reaching the Earth 
change the energy input that forces atmospheric composition 
and temperature and the Earth’s climate (for a recent review, 
see Gray et al., 2010). The spectral distribution of solar irra- 
diance is vital for determining the Earth’s response to solar I45 
cycle variations (Zhong et al., 2008). Although the 1 1 -year 
solar cycle relative variation in TSI is only 0.1%, the solar 
variation of the SSI, as a function of wavelength, can be 
much larger — greater than 10% in the oxygen Schumann- 
Runge bands, for example. The absolute magnitude of the , 60 
solar radiation in the UV is much smaller than that at longer 
wavelengths, but its absorption in the middle atmosphere is 
critical to maintaining the composition, thermal structure, 
and dynamics of the atmosphere. 

With the advent of space-based observations, beginning , 55 
in 1978, the variability of the solar irradiance on minute-to- 
decadal timescales has been measured (Frohlich and Lean, 
2004). The traditional view of the 1 1-year solar cycle holds 
that the UV, visible, and IR portions of the SSI rise and fall in 
phase with the TSI and other solar proxies, such as the solar , 6 o 
10.7-cm radio flux, F 10 .7 (Lean, 2000). This view is based 
on both observations and our understanding of solar physics. 
How changes in the solar energy reaching the Earth impact 
climate (temperature and precipitation) variability has been 
the focus of intense investigation (e.g., Haigh, 1996; Shindell t6 5 
et al, 1999, 2001 ; Rind, 2002; Rind et al., 2004, 2008; Ineson 
et al., 2011). 

Accurate characterization of the SSI is so important that 
not one but two separate instruments measuring UV SSI 
were flown on the Upper Atmosphere Research Satellite 170 
(UARS): the SOLar STellar Irradiance Comparison Exper- 
iment (SOLSTICE) (Rottman et al„ 1993) and the Solar Ul- 
traviolet Spectral Irradiance Monitor (SUSIM) (Brueckner 
et al., 1993). Recent satellite measurements by SORCE, 
however, challenge the canonical view of SSI solar cycle 175 
variability represented by the work of Lean (2000). Since 


it was launched in 2003, the SORCE mission has been mea- 
suring SSI and TSI (Rottman et al., 2005; McClintock et al., 
2005; Kopp et al., 2005; Kopp and Lean, 201 1 ). 

Solar UV radiation has two primary direct effects on 
the atmosphere: atmospheric heating and the photolysis of 
molecular oxygen (O 2 ) and trace gases, such as ozone (e.g., 
Brasseur and Solomon, 1986). In the stratosphere and meso- 
sphere, solar heating is principally through absorption of UV 
radiation by ozone in the Hartley (200-300 nm) and Hug- 
gins (300-360 nm) bands, with absorption of visible light by 
ozone becoming more important in the lower stratosphere. 
Absorption by the 0 2 Schumann-Runge system (<205 nm) 
is important in the upper mesosphere, but heating via the 
Herzberg continuum (190-242 nm) has only a second-order 
effect in the middle atmosphere. Absorption of visible and 
near-infrared radiation in the lower atmosphere by ozone, 
ILO, CO 2 , NO 2 , and other minor gases contributes to heat- 
ing in this region. Through photolysis, solar radiation is the 
driver of atmospheric chemistry; in the stratosphere, solar ra- 
diation is directly responsible for the ozone production and 
loss processes that together control the amount of ultraviolet 
radiation reaching the lower atmosphere and surface. 

Precisely how ozone and temperature respond to solar cy- 
cle variability depends on the spectral dependence of SSI 
variation and the interplay between atmospheric heating and 
photolysis. Direct heating increases during the ascending 
phase of the solar cycle, increasing atmospheric temper- 
ature and decreasing ozone concentrations, due to faster 
temperature-dependent ozone loss reactions, which conse- 
quently cause a negative feedback on temperature as the pri- 
mary absorber of radiation is destroyed. Increased solar ra- 
diation also increases ozone production (through O .2 photol- 
ysis) and loss (through increased atomic oxygen concentra- 
tions resulting from greater ozone photolysis). Net ozone 
production leads to more absoiption and heating. The bal- 
ance of these interrelated and competing processes and their 
wavelength dependencies can be explored using models. 

Typical atmospheric general circulations models (GCMs) 
include the solar direct atmospheric heating term and have 
interactive radiation and dynamics, but chemistry (including 
photochemistry) is specified and does not include solar mod- 
ulation of chemistry or chemistry feedbacks through changes 
in composition. Chemistry-transport models (CTMs) in- 
clude changes in photochemistry, but the radiation and dy- 
namics are specified and non-interactive with the chem- 
istry. What is therefore needed to address the effects of 
changing solar irradiance on both heating and photolysis is 
a chemistry-climate model, which treats both effects self- 
consistently, with interactive chemistry, radiation, and dy- 
namics. 

Several modeling studies have examined the impact of the 
solar cycle on the middle atmosphere, using 2-D models 
(e.g., Garcia et al., 1984; Fleming et al., 1995; Huang and 
Brasseur, 1993; Haigh et al., 2010), 3-D CTMs (e.g., Sto- 
larski et al., 2006), and 3-D CCMs (e.g.. Marsh et al., 2007; 
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Austin et al., 2008). These previous modeling studies have 
been based at least in part on the NRL reconstruction of SSI 
and have found temperature and ozone responses to the 1 1-230 
iso year solar cycle that, although different from measurement- 
derived responses in certain details, are largely consistent 
with long-term observations (e.g., Austin et al., 2008). The 
ozone response, for example, has been found to be in-phase 
with TSI in the stratosphere and lower mesosphere. 23 s 

is5 More recently, studies have looked at the SORCE SSI 
specifically. Cahalan et al. (2010) investigated the impact 
of the new SORCE SSI in the Goddard Institute for Space 
Science (GISS) modelE GCM in comparison to the NRL re- 
construction of SSI and found that the SORCE SSI led to 240 
iso stratospheric heating that agreed better with the solar cy- 
cle impact on temperature inferred from Halogen Occulta- 
tion Experiment (HALOE) measurements (Remsberg, 2008). 
Haigh et al. (2010) modeled both the NRL and SORCE SSI 
in a 2-D model and predicted ozone variations in phase with 245 
las the solar cycle below 45 km and out of phase above. They 
also included observational evidence to support the SORCE 
SSI fingerprint in stratospheric ozone using Aura/Microwave 
Limb Sounder (MLS) observations. The period of avail- 
able SORCE data and the model comparison of Haigh et al. 2 S 0 
200 with MLS covered only part of a solar cycle, however, mak- 
ing it very difficult to draw statistically robust conclusions 
(Garcia, 2010). And although both studies led to conclu- 
sions supporting the SORCE SSI, they predicted opposite 
surface temperature responses, with Cahalan et al. in phase 
see with TSI and Haigh et al. out of phase with TSI. Merkel 
et al. (2011) performed simulations using the Whole At- 
mosphere Community Climate Model (WACCM), predict- 
ing an ozone response similar to Haigh et al. when using 
SORCE SSI, out-of-phase with the solar cycle above ap- 
210 proximately 40 km. They performed a regression analysis of 
Thermosphere Ionosphere Mesosphere Energetics and Dy- 
namics (TIMED)/Sounding of the Atmosphere using Broad- 
band Emission Radiometry (SABER) ozone data, a record 
somewhat longer (8 years) than the MLS time series analyzed 
215 by Haigh et al., and also inferred an out-of-phase response in 
the lower mesosphere. 

265 

How the atmosphere responds to solar cycle variations in 
the spectral distribution of SSI has important implications 
for the projection of atmospheric composition and climate 
220 on decadal and longer timescales, and even independent of 
the current debate regarding the magnitude and sign of the ^ 
solar cycle in SSI, a better understanding of the atmospheric 
response to changes in SSI is needed. 

3 Solar spectral irradiance 
225 3.1 NRL SSI 

275 

Realistic simulations of the solar-forced climate on multi- 
decadal timescales rely on long-term solar irradiance recon- 


structions. The NRL SSI reconstruction used in this work 
(Lean, 2000) was derived using a parameterization of the 
observed temporal variation of SSI over the past 30 years 
with proxy indicators needed to represent the known wave- 
length sources of variability in SSL Historical records of 
these proxies were then used to infer SSI variations before 
satellite observations were available. The model of total so- 
lar irradiance variations used was described by Frohlich and 
Lean (2004). The specific parameterizalions used were those 
of Frohlich and Lean ( 1 998) and Lean (2000), which used 
the Physikalisch-Meteorologisches Observatorium Davos 
(PMOD) multi-satellite composite. This self-consistent re- 
construction of TSI and SSI includes both annual and 
monthly total solar irradiance and spectral solar irradiance 
from near hydrogen Lyman-o (120nm) through 100 pm, 
from January 1882 through December 2008, comprising 12 
solar cycles. The irradiance values have been scaled for the 
present analysis by 0.9965 to match the SORCE/Total Irradi- 
ance Monitor (TIM) absolute scale, which has been verified 
by NIST to be more accurate (see also Kopp and Lean, 2011). 
Note that the NRL SSI reconstruction is widely used in cli- 
mate simulations, including those in the IPCC reports (e.g., 
IPCC, 2007). 

3.2 SORCE SSI 

The SORCE mission comprises the SOLar STellar Ir- 
radiance Comparison Experiment (SOLSTICE), the 
Spectral Irradiance Monitor (SIM), TIM, and the XUV 
Photometer System (XPS) (Rottman, 2005). Publicly 
available SOLSTICE data, at wavelengths shorter than 
310nm, were used in this work. SIM measures SSI 
from 200 to 2270 nm, although only data at wavelengths 
longer than 310nm are currently publicly available (see 
http://lasp.colorado.edu/sorce/data/data-product_summary.htm). 
The SORCE team, however, provided us with SIM spectral 
time series data from 200 to 1600nm (J. W. Harder, per- 
sonal communication) along with analysis procedures to 
remove data spikes resulting from spacecraft roll maneuvers. 
Following the recommendation of the SORCE team, SOL- 
STICE data were used at wavelengths less than 244 nm, as 
the signal-to-noise ratio of the SIM measurements is poor in 
this wavelength region. SORCE data are available starting 
in April 2004, in the middle of the declining phase of solar 
cycle 23, extending through the most recent solar minimum, 
near the end of 2008, and to present. 

3.3 Derivation of solar maximum/minimum cases 

All the model simulations presented in this paper were per- 
formed in time-slice or non-time-dependent mode, stationary 
in time with respect to solar forcing. Therefore we consid- 
ered two pairs of SSI spectra representative of solar maxi- 
mum (S R , ax) and minimum (S m j„) conditions, derived from 
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the NRL and SORCE SSI datasets. The spectral differences 
are shown in Fig. 1 . 330 

For the NRL SSI cases, we chose monthly mean spec- 
tra from the reconstruction representative of S max (August 
2001) and S min (June 2007) conditions. These two partic- 
ular months were selected because they were close to solar 
max/min conditions, lie on a linear fit correlating the NRL 335 
SSI in the ultraviolet and F 10 . 7 , and are separated by nearly 
100 units of F 10 . 7 . This allows the results to be interpolated 
or extrapolated approximately to other solar activity condi- 
tions, as the effects scale linearly over this range of magni- 
tudes (not shown). 340 

Because SORCE SSI measurements cover only a part of 
the spectral range simulated in the models used in this work, 
the following procedure was used to develop S ma x(SORCE) 
and S m j n (SORCE) SSI spectra. As mentioned above, SOL- 
STICE data were used at wavelengths below 244 nm and 345 
SIM was used from 244 to 1600nm (SIM actually mea- 
sures the SSI to roughly 2270 nm, but the inferred solar cy- 
cle is considered to be more reliable below 1600nm). At 
longer wavelengths, the spectra were supplemented with the 
NRL SSI, but middle atmosphere ozone and temperature 350 
are insensitive to these longer wavelengths. This compos- 
ite SSI time series was correlated with Fj 0.7 from Septem- 
ber 2004 through November 2007 on a per-wavelength ba- 
sis, fx(Fio,j). S nl j n (SORCE) conditions were then con- 
structed from the fit(s) at solar minimum conditions (Novem- 355 
ber 2007); S nlllx (SORCE) SSI was inferred for solar condi- 
tions of 100 units F 10.7 greater as /A(Fio. 7 (S m i n + 100 )). This 
linear fitting approach allowed us to infer solar maximum 
conditions based on the available SORCE time series. 


4 Chemistry-climate models 
4.1 GEOS CCM 

365 

The 3-D model used in this study was the Goddard Earth 
Observing System chemistry-climate model (Pawson et al., 
2008, and references therein). The CCM uses the GEOS-5 
atmospheric general circulation model (Rienecker et al., 
2008), version 5.2.0, coupled with the StratChem chemical 
solver developed as a part of the GSFC 3-D CTM (Douglass 
et al., 1996; Pawson et al., 2008). Eyring et al. (2006) showed 370 
that GEOS CCM simulations agree well with observations, 
both in terms of dynamics/transport and chemistry, and that 
the evaluated quantities lie well within the range of results 
obtained by all the models in the study (see also Austin et al., 

2008). 375 

The GEOS CCM traditionally uses a fixed input solar 
spectrum, representative of mean solar cycle conditions, and 
has in fact been used as a no-solar cycle reference model 
in past CCM intercomparisons (Austin et al., 2008). For this 
work, it was therefore necessary to implement the solar cycle aso 
in the model. In the GEOS CCM, the solar forcing of atmo- 


spheric heating and photolysis are treated separately. This 
structural aspect of the model code was beneficial for the 
present study, as we wanted to examine solar cycle effects 
on heating and photolysis independently. 

Six time-slice simulations were performed using the 
GEOS CCM forced with the NRL SSI reconstruction and 
solar cycle flux conditions listed in Table 1 (first six table 
entries) in order to evaluate the independent and combined 
effects of SSI variations through photolysis and direct atmo- 
spheric heating. We replaced the default solar spectrum with 
S max and S nl in conditions from the NRL SSI. Two additional 
time-slice simulations were performed based on the SSI mea- 
sured by SORCE. 

Each time-slice run was a 30-year simulation in which 
greenhouse gases (other than CO 2 ) and ozone-depleting sub- 
stance emissions were fixed at 2005 values by specifying 
their mixing ratios in the lowest model layer, based on es- 
tablished inventories (World Meteorological Organization, 
2003; Eyring et al., 2006). Aerosols, including sulfate sur- 
face area, were fixed to 1979 values. In contrast, sea surface 
temperatures (SSTs) and sea ice distributions were allowed 
to follow their observed 1979-2008 evolutions, specified in 
the model using the Hadley Centre Ice and Sea- Surface Tem- 
perature dataset (Rayner et al., 2003). For consistency with 
the climate change signal in SSTs, CO 2 was also allowed to 
increase over the same time period according to the IPCC 
A1B scenario (IPCC, 2000), with the knowledge that SSTs 
increase with the atmospheric CO 2 burden. All the model 
results shown in this work are in the form of S max — S m j n rel- 
ative differences, so changing CO 2 and SSTs in otherwise 
time-slice simulations would have only a second-order effect 
and should not bias the present results. All CCM simula- 
tions started in SST/C0 2 month January 1979 from a pre- 
viously spun-up CCM run using mean solar flux conditions. 
For all the comparisons shown here, the first five model years 
have been discarded, which would exclude the period of at- 
mospheric adjustment to S max /S m i n conditions. The response 
in total ozone indicated that the relaxation time of the model 
atmosphere was actually of the order of months. 

4.2 GSFC 2-D chemistry-radiation-dynamics coupled 
model 

We also performed a series of simulations using the God- 
dard Space Flight Center 2-D coupled chemistry-radiation- 
dynamics model (see Table 1 ). With its computational effi- 
ciency, the GSFC 2-D model is ideal for studying the sensi- 
tivity of the atmosphere to multiple perturbations in the input 
SSI and is better suited for conducting a large number of ex- 
periments. Two-dimensional models, including the GSFC 2- 
D model, have been used extensively in international assess- 
ments of the stratosphere (e.g., World Meteorological Orga- 
nization, 2011), and the GSFC 2-D model has been used pre- 
viously to study the atmospheric response to both 27-day and 
1 1-year solar cycle variability (Fleming et al., 1995, 2007). 
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Because its chemistry, radiation, and dynamics are fully 
coupled, the GSFC model is essentially a 2-D CCM. The 
current version of the GSFC 2-D was described by Fleming 
385 et al. (2011), who have shown that the model simulations 435 
compare very well with stratospheric observations and the 
GEOS 3-D CCM. For consistency with GEOS CCM simula- 
tions presented here, the same input SSI spectra representing 
S max and S m in conditions derived from the NRL and SORCE 
390 SSI were used. Thirty-year time-slice runs were performed «o 
with the 2-D model, and the final year of simulation was an- 
alyzed (averaging over more years is not necessary because 
of the almost negligible interannual variability in the 2-D 
model, compared with the GEOS CCM). 

445 

396 5 Discussion 

5.1 Solar cycle response through heating and photolysis 

The atmosphere responds to the 1 1-year cycle in solar radia- 460 
tion through direct atmospheric heating and the photolysis of 
0 2 and trace gases, such as ozone. These different manifolds 
400 have been examined in a GCM by performing simulations 
alternatively with the solar cycle in direct radiative heating 
and with the solar signal in specified ozone variations, re- 
spectively (Shibata and Kodera, 2005). The advantage of us- 
ing a CCM, however, is that secondary feedback mechanisms 455 
4os are coupled. For example, when the solar cycle is included 
only in heating (e.g., runs S max H and S^H in Table 1), al- 
though the solar cycle in photolysis is intentionally excluded, 
the ozone loss due to greater direct heating at S max and its 
feedbacks are still captured. 460 

410 The response of the atmosphere through heating and pho- 
tolysis were examined individually and in combination in the 
GEOS CCM by including the NRL SSI solar cycle selec- 
tively in direct heating and photolysis (and both combined), 
with the model runs listed in Table 1 . All results are reported 465 
415 relative to 100 units of F10.7 (the total F10.7 variation in a 
typical solar cycle is about 125). 

5.1.1 Ozone 

470 

The solar cycle impact on the ozone profile is considered 
first. Fig. 2 shows profiles of the Smax— Smin ozone relative 
420 difference, averaged over 60°S-60°N, based on 25 years of 
time-slice simulation. In the heating+photolysis simulation, 
SmaxPH— S m i n PH, the magnitude of the solar cycle difference 
peaks in the upper middle stratosphere near 3 hPa at about 
+2.0% per 100 units of F10.7. This is generally consistent 475 
425 with other modeling efforts and observations (e.g., Austin 
et al., 2008). The heating-only case shows a slight negative 
ozone response above 30hPa, due to increased ozone loss 
rates at warmer temperatures, with a positive response in the 
lower stratosphere where dynamics plays a larger role. The4go 
430 photolysis-only simulation captures almost all of the solar 
cycle effect on ozone. Comparison of the three simulations 


also reveals that the direct heating and photolysis response 
mechanisms are roughly linearly separable above 30 hPa. 

The solar cycle in SSI leads to a mean total ozone en- 
hancement of 3 Dobson Units per 100 units of F\ 0.7 in the 
tropics going from S m j„ to S max (see Fig. 6). This overall 
effect is similar to previous modeling efforts (e.g., Austin 
et al., 2008; Marsh et al.. 2007) and the solar cycle term de- 
rived from regression analyses of merged ozone data (e.g., 
Stolarski et al., 2006) and ground-based measurements (e.g., 
Randel and Wu, 2007). At polar latitudes, the modeled inter- 
annual variability is quite significant, with larger variability 
in the Northern Hemisphere (not shown). Similar conclu- 
sions can be drawn from the total ozone response regarding 
the relative contributions of direct heating and photolysis and 
their linear separability (not shown). 

Thus, in the case of ozone, the solar cycle response is 
caused almost entirely through photolysis, with heating pro- 
viding only a small, mostly negative effect. This is impor- 
tant in that it suggests that CTMs (e.g., Stolarski et al., 2006) 
should be sufficient to capture the response in ozone to so- 
lar cycle variability — chemistry feedback to the radiation and 
dynamics is of secondary importance. 

5.1.2 Temperature 

The corresponding plot for the solar cycle impact on tem- 
perature is shown in Fig. 3. The maximum effect occurs 
near the stratopause at about 0.6 K per 100 units of F10.7, 
at 6Q°S-60°N. In contrast to ozone’s response to the solar 
cycle, however, both direct heating and photolysis contribute 
positively to warmer temperatures at solar maximum. In fact, 
the additional ozone produced via O2 photolysis with greater 
UV flux at S max leads to a larger temperature response than 
the solar cycle in direct heating itself. Further, the two mech- 
anisms are close to linearly additive throughout most of the 
stratosphere and lower mesosphere. 

The finding that direct heating and photolysis both make 
significant contributions to the increased stratospheric tem- 
peratures associated with solar maximum conditions implies 
that GCMs do not adequately represent the atmospheric re- 
sponse to SSI solar cycle variability. Without including the 
solar cycle in photolysis, a GCM would substantially under- 
estimate the temperature response. 

5.2 Response to the SORCE SSI solar cycle 

Using the S max , S m j n SSI cases described in Sect. 3.3, we next 
compare the atmospheric responses to the NRL and SORCE 
SSI (see SSI spectra in Fig. I ). The results for ozone are 
shown in Figs. 4 and 5. The response using the NRL SSI 
shown in Fig. 4 is the same as the heating+photolysis case 
in Fig. 2. The modeled response shown in Fig. 5 is statisti- 
cally significant (95% confidence level) throughout most of 
the plotted range, as indicated. The ozone response to the 
NRL SSI is positive at all altitudes, up through the middle 
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mesosphere. The response using the SORCE SSI is very dif- 
ferent in comparison. Although the ozone response in most 
485 of the stratosphere going from S m j n to S max is positive, ors 4 o 
in phase with the solar cycle in TSI, above about 2 hPa the 
response is actually negative (up to 0.2 hPa). This modeled 
response is largely consistent with previously published re- 
sults (e.g,, Haigh et al., 2010; Merkel et al., 2011) and with 
4 so the GSFC 2-D model (see Fig. 9). 545 

Recent analyses of ozone trends over the SORCE time 
period, representing part of solar cycle 23, have been used 
to infer an out-of-phase ozone response in the upper strato- 
sphere/lower mesosphere from Aura/MLS (Haigh et al., 

495 2010) and TIMED/SABER (Merkel et al., 2011) measure -550 
ments. It must be stressed, however, that extracting statis- 
tically significant ozone trends from partial solar cycles is 
extremely challenging (Garcia, 2010). The ozone signal is 
small, and there are other competing variations on various 
500 time scales, such as those caused by the quasi-biennial oscil- 
lation, the El Nino Southern Oscillation, volcanic eruptions, 
declining concentrations of ozone-depleting substances, and see 
climate change (e.g., Stolarski et al., 2006). 

Longer satellite ozone records are available and have been 
505 studied as well, in detail. The solar cycle in ozone inferred 
from three solar cycle’s worth of satellite observations by the 
Solar Backscatter Ultraviolet Instrument (SBUV), the Strato- sm 
spheric Aerosol and Gas Experiment (SAGE), and HALOE, 
also shown in Fig. 4, lead to a combined ozone response that 
510 is positive everywhere in this altitude range (Austin et al., 
2008; Soukharev and Hood, 2006), similar to the NRL SSI- 
based simulation. These longer-term studies, however, do 555 
not include the declining phase of solar cycle 23, and thus a 
unique solar cycle variation during the SORCE mission (e.g., 

515 2004-present) cannot be ruled out. Another complication in 
these long satellite records is that occultation instruments, 
such as SAGE and HALOE, make their observations near the 570 
terminator, at a time of rapid diurnal variation in mesospheric 
ozone (e.g., Merkel et al., 201 1 ). 

62 o The longest, perhaps least ambiguous record of the ozone 
response to the solar cycle is that of total column ozone, 
as shown in Fig. 6. The solar cycle response derived from 575 
ground-based observations (Randel and Wu, 2007) agrees 
well with the GEOS CCM using NRL SSI. The inferred so- 
525 lar cycle from merged Total Ozone Mapping Spectrometer 
(TOMS) and SBUV total ozone data (Stolarski et al., 2006; 
Randel and Wu, 2007) is larger than that from ground-based seo 
measurements (roughly double in the tropics) but is similar to 
the modeled response based on SORCE SSI. Unfortunately, 

530 large statistical uncertainties make it difficult to distinguish 
between these four representations of the solar cycle in total 
ozone. 565 

The modeled solar cycle response in temperature is shown 
in Figs. 7 and 8. In the NRL SSI case, the temperature re- 
535 sponse to the solar cycle is roughly 0.6 K per 100 units of 
Fm.r in the tropics, peaking near the stratopause. Note that 
this response in the GEOS CCM, which includes the effect 590 


of feedback from net solar cycle ozone production via O 2 
photolysis (and thus additional heating) is significantly larger 
than that predicted with a GCM without coupled chemistry 
(e.g., Cahalan et al., 2010). The modeled response using the 
SORCE SSI is much larger still, almost 2 K per 100 units of 
f\o. 7 , peaking near 2 hPa. The fact that the response peaks 
at a lower altitude is not surprising, considering the relative 
ozone loss in the lower mesosphere and the large SORCE 
SSI variation in the ozone Huggins bands wavelength re- 
gion (300-360 nm), which penetrates lower into the atmo- 
sphere. At its peak, the SORCE SSI case is over twice as 
large as the solar cycle temperature response inferred from 
either HALOE (Remsberg, 2008) or microwave temperature 
sounding (Fig. 7; Austin et al., 2008). 

5.3 Decomposition of the ozone response as a function 

of wavelength 

The response of the atmosphere to the 1 1-year solar cycle is 
sensitive to the spectral characteristics of the SSI solar cy- 
cle variation. As discussed in Sect. 2, increasing UV radi- 
ation from S min to S max increases both O 2 photolysis, lead- 
ing directly to ozone production, and ozone photolysis to 
form atomic oxygen, which leads to catalytic ozone loss. 
Although ozone photolysis occurs between roughly 200 and 
340 nm, the ozone balance is dominated by O 2 photolysis 
below 242 nm. At wavelengths longer than this O 2 photoly- 
sis threshold, the situation is driven by ozone photolysis and 
catalytic loss. This means that the net ozone production/loss 
in the photochemically controlled upper stratosphere/lower 
mesosphere is determined by the relative contributions of 
these two wavelength ranges. Coincidentally, the 242-nm 
threshold is close to the recommended minimum wavelength 
for use of the SIM SSI data (see SSI spectra in Fig. 1 ), so the 
net ozone response to the solar cycle results from the relative 
size of the solar cycle inferred from the SOLSTICE and SIM 
instruments of the SORCE mission. 

To study the details of the contributions of the SSI varia- 
tions on either side of the 242-nm threshold to the ozone re- 
sponse, we used the GSFC 2-D model, described in Sect. 4.2. 
The 2-D model, with fully coupled chemistry-radiation- 
dynamics, is essentially a 2-D CCM. Although the 2-D 
model may not capture all the details of ozone transport in 
the lower stratosphere, it is more than adequate for the pho- 
tochemical problem posed here. Fig. 9 shows the ozone 
response to the solar cycle as described by the NRL and 
SORCE SSI and the excellent agreement between the 2-D 
model and the GEOS CCM. This figure also shows 2-D 
model simulations where the solar cycle in SSI is applied al- 
ternatively at wavelengths less and greater than 242 nm. The 
ozone production-dominant short-wavelength region of the 
solar cycle leads to enhanced ozone throughout the strato- 
sphere and lower mesosphere, whereas the solar cycle at 
longer wavelengths leads to ozone loss at virtually all alti- 
tudes. 
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The ozone responses at wavelengths less and greater 
than 242 nm are qualitatively similar for both the NRL and 
SORCE SSI. The quantitative difference in the total re- 646 
sponses between the two SSI specifications results from the 
595 balance of the two wavelength regimes. The magnitude of 
the solar cycle in the SOLSTICE part of the SSI spectrum 
(A <242 nm) is greater than the NRL SSI, leading to a peak 
ozone production more than a factor of 2 larger. The mag- eso 
nitude in the SIM region (A >242 nm), however, is greater 
coo by roughly a factor of 6 at its peak. The ozone response 
to the solar cycle in both parts of the spectrum combined is 
almost exactly equivalent to the linear combination of the re- 
sponses in the individual wavelength ranges (not shown ex- ess 
plicitly). And it is the balance of these ozone production and 
60 s loss regimes that results in the very different total responses 
in the model when forced with the NRL and SORCE SSI, 
including the altitude at which the SORCE ozone response 
goes from in-phase (with respect to the solar cycle in TSI) in eeo 
the lower and middle stratosphere to out-of-phase near and 
eio above the slratopause. 

6 Conclusions 665 

Model simulations using the GEOS CCM and GSFC 2-D 
coupled model have been used to further clarify how the at- 
mosphere responds to the SSI solar cycle. Using the GEOS 
sis CCM to explore the effects of SSI variations on direct at- 
mospheric heating and photolysis both separately and when 670 
combined, we have quantified their relative contributions in 
relaying solar cycle variations to stratospheric ozone and 
temperature. We have shown how the response of ozone to 
62 g the solar cycle is almost entirely due to the solar cycle in 676 
photolysis, whereas the solar cycle in both photolysis and di- 
rect atmospheric heating are important for the temperature 
response. This means that CTMs should simulate the ozone 
response well, while GCMs without coupled chemistry will 
ess underestimate the temperature response significantly. Also, 
the responses in both ozone and temperature from photolysis 
and heating are linearly independent to first order. 

Simulations presented here with the GEOS CCM and 
GSFC 2-D model confirm recent simulations based on both 
630 the NRL and SORCE SSI (Haigh et al., 2010; Merkel et al., 
2011). The NRL SSI produces an ozone response that is in 
phase with the solar cycle in TSI throughout the stratosphere ess 
and lower mesosphere, whereas the SORCE SSI produces 
an out-of-phase response in the uppermost stratosphere and 
635 lower mesosphere. Although this model result based on the 
SORCE SSI is more consistent with the ozone solar cycle in- 
ferred from recent satellite observations (Haigh et al., 2010; 
Merkel et al., 2011), caution should be exercised when draw- 
ing conclusions. Using climate or chemistry-climate model 
640 output in comparison with a solar cycle inferred from a short 
data record is probably not the best way to validate the solar 69f , 
cycle inferred from SSI measurements. Even the most com- 


prehensive models are subject to numerous simplifications. 
In addition, inferring solar cycle variations from part of a so- 
lar cycle is challenging (Garcia, 2010). A more convincing 
approach is to make comparisons with the SSI observations 
themselves, such as has been done by Morrill et al. (2012) 
and DeLand and Cebula (2012). Ultimately, continued ob- 
servation of the SSI by SORCE and subsequent missions will 
likely be necessary to definitively resolve the NRL-SORCE 
discrepancies. Further, it is possible that interpretation of the 
SIM measurements will change as more of solar cycle 24 is 
observed. 

We have also shown using 2-D model simulations that the 
net ozone response in the middle atmosphere results from 
the SSI-determined balance of ozone production at wave- 
lengths less than 242 nm and ozone destruction at longer 
wavelengths. These linearly separable wavelength regimes 
roughly correspond to those of the SOLSTICE and SIM 
instruments on SORCE, respectively, meaning the relative 
sizes of the UV SSI solar cycle variations inferred from these 
two instruments determine the magnitude and phase (with re- 
spect to the solar cycle in total solar irradiance) of the ozone 
response. An analysis of the atmospheric response at finer 
wavelength resolution and its linearity could allow for a bet- 
ter prediction of the atmospheric response to arbitrary SSI 
variations. 
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Fig. 1. SSI spectral differences, S max relative to S m ,n per 100 units of 
Fio. 7, derived from NRL and SORCE SSI time series. The vertical 
dashed line at 244 nm represents both the wavelength cut-off sep- 
arating the SORCE SOLSTICE and SIM spectra used and the O2 
photolysis threshold (O2 photolysis occurs at wavelengths shorter 
than 242 nm). 



Ozone Relotive Difference (%/100 units F 107 ) 


Fig. 2. Ozone response to the 1 1-year solar cycle, Sm ax — S m i„ (rela- 
tive %) per 100 units of F10.7, in the stratosphere and lower meso- 
sphere and averaged over 60°S-60°N latitude, as modeled by the 
GEOS CCM. Three cases are shown, where the solar cycle has been 
applied (1) only to direct radiative heating, (2) only to photolysis, 
and (3) with both effects combined. 
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Table 1. List of CCM time-slice simulations performed with the GEOS CCM and GSFC 2-D model 


Designation 

Model 

SSI Source 

Solar Conditions 

Solar Cycle in 

SmaxH 

GEOS CCM 

NRL 

Max 

Heating 

SminH 

GEOS CCM 

NRL 

Min 

Heating 

SmaxP 

GEOS CCM 

NRL 

Max 

Photolysis 

S„„„P 

GEOS CCM 

NRL 

Min 

Photolysis 

Smax PH 

GEOS CCM 

NRL 

Max 

Heating+Photolysis 

SmmPH 

GEOS CCM 

NRL 

Min 

Heati ng+Photoly sis 

SmaxPH(SORCE) 

GEOS CCM 

SORCE 

Max 

Heating+Photolysis 

S m ,„PH(SORCE) 

GEOS CCM 

SORCE 

Min 

Heating+Photolysis 

2DS raas PH 

GSFC 2-D 

NRL 

Max 1 

Heating+Photolysis 

2DSminPH 

GSFC 2-D 

NRL 

Min 

Heating+Photolysis 

2DS mas PH(SORCE) 

GSFC 2-D 

SORCE 

Max 1 

Heating+Photolysis 

2DS mi „PH(SORCE) 

GSFC 2-D 

SORCE 

Min 

Heating+Photolysis 


1 S max conditions in the GSFC 2-D model were applied alternatively at wavelengths <242 nm, >242 nm, and over the entire spectrum (see 
Fig. 9). 



Ozone Relotive Difference (%/100 units F t07 ) 



Temperature Difference (K/100 units F I07 ) 


Fig. 3. Same as Fig. 2, but for the temperature response to the solar 

cycle. 


Fig. 4. Simulated ozone solar cycle response in the GEOS CCM 
using both NRL and SORCE SSI, averaged over 25°S-25°N lati- 
tude. Also shown is the solar cycle variation in ozone deduced from 
SBUV, SAGE, and HALOE satellite measurements (Austin et al. 
(2008), Fig. 4; see also Soukharev and Hood (2006)). 




W. H. Swartz et al.: Middle atmosphere response to the solar cycle in SSI 


(a) NRL SSI (b) SORCE SSI 



60S 30S EQ 30N 60N 60S 30S EQ 30N 60N 

Latitude (’) Lotitude (*) 


Fig. 5. Simulated ozone solar cycle response in the GEOS CCM, similar to Fig. 4 but in altitude-latitude cross section, using (a) NRL and 
(b) SORCE SSI. The colored (non-white) portions of the plots indicate regions where the ozone response is statistically significant at the 
95% confidence level. 



60S 30S EQ 30N 60N 

Latitude (*) 


Fig. 6. Total ozone response to the 1 1-year solar cycle, as modeled 
by the GEOS CCM and inferred from observations by regression 
analyses of the TOMS. SBUV merged total ozone dataset (Stolarski 
et al., 2006, Fig. 6) and ground-based observations (Randel and Wu, 
2007, Fig. 12). 



Temperature Difference (K/100 units F t07 ) 


Fig. 7. Same as Fig. 4, but for temperature. SSU=Stratospheric 
Sounding Unit; MSU=Microwave Sounding Unit; data from 
(Austin et al., 2008, Fig. 12). 
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(a) NRL SSI 



Latitude (*) 



60S 30S EO 30N 60N 

Latitude (*) 


Fig. 8. Same as Fig. 5 but for temperature. 


(a) NRL SSI 



Relative Difference (%/100 units F 107 ) 


(b) SORCE SSI 



Relative Difference (%/100 units F 107 ) 


Fig. 9. Simulated ozone response to the solar cycle at wavelengths less and greater than 242 nm using both (a) NRL and (b) SORCE SSI in 
the GSFC 2-D model, along with the total responses, averaged over 60 C S-60°N latitude. Also shown are the analogous GEOS CCM total 
responses (corresponding to the simulations in Fig. 4 but averaged over 60°S-60°N) for comparison. 
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